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Erosion  Damage  in  Carbon-Carbon  Composites 
at  Hypersonic  Impact  Velocities 


William  F.  Adler* 
and 

*■. 

Anthony  G.  Evans 


Abstract 


The  damage  to  three-dimensional  orthogonal  carbon-carbon  composites 
created  by  particle  Impacts  ranging  from  1800  to  6000  ms~^  has  been  char- 
acterized for  three  weave  constructions  at  both  ambient  and  at  temperatures 
up  to  3600°C.  The  subsurface  damage  modes  In  the  vicinity  of  the  Impact 
crater  are  described  for  the  range  of  conditions  Investigated.  Hypotheses 
are  advanced  concerning  the  failure  modes  and  kinetics  for  the  material 
ejected  from  the  crater  based  on  examination  of  sectioned  specimens  and 
high-speed  photographic  sequences.  Penetration,  compliance,  and  plane- 
strain  compression  tests  are  utilized  to  supply  mechanical  characteri- 
zation data  for  the  dominant  modes  of  failure, which  are  analyzed  In  terms 
of  the  dynamic  stress  states  which  develop  In  the  three-dimensional 
composites  during  the  collision  cycle. 


*Dlrector,  Materials  Science  Section,  Effects  Technology,  Inc.,  5383 
Hollister  Ave.,  Santa  Barbara,  CA  93111. 

Group  Leader,  Structural  Ceramics , Structural  Materials  Department, 
Rockwell  International  Science  Center,  1049  Camlno  Dos  Rios,  Thousand 
Oaks,  CA  91360. 
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Investigation  of  the  erosive  response  of  fibrous  composites  exposed 
to  natural  (raln/lce)  and  Induced  (dust)  particulate  environments  Is  stimu- 
lated by  the  need  for  re-entry  vehicle  protection  materials.  The  erosion 
performance  of  existing  system  materials  has  been  evaluated  using  ground 
Intact  simulation  testing  and  verification  flight  tests. In  the  eval- 
uation of  materials  fabricated  for  Improved  erosion  resistance,  all  material 
evaluation  has  been  directly  related  to  erosion  test  performance.  Little 
attention  has  been  given  to  Information  for  material  development  activities. 
It  has  thus  been  difficult  for  material  manufacturers  to  relate  erosion 
test  results,  such  as  mass  loss  ratios,  to  material  construction  or  pro- 
cessing variables. 


Mass  loss  determinations  have  been  performed  either  by  single  particle 

Impact  procedures or  by  continuous  erosion  testing  In  ballistic  range 

(1-3) 

facilities.  Mass  loss  laws  have  then  been  derived  for  specific  materials 

based  on  experimental  correlations.  However,  the  single  particle  studies 
have  not  adequately  represented  effects  of  multiple  collisions,  while  the 
laws  derived  from  the  continuous  erosion  studies  have  focused  on  the  re- 
sult and  not  the  process  of  erosion.  In  both  cases  extrapolation  to  service 
conditions  has  been  made  from  a limited  data  base  and  without  consideration 
of  the  physical  processes  by  which  cratering  occurs.  Recently,  attempts 
have  been  made  to  Interrelate  all  aspects  of  the  erosion  process, but 
the  extensive  use  of  empirical  variables  has  limited  the  utility  of  this 
procedure  for  the  a priori  prediction  of  erosive  material  removal. 


Each  single  particle  Impact  causes  substantial  subsurface  damage 
In  addition  to  mass  loss.^^^^  Furthermore,  after  many  sequential  single 
particle  Impacts  on  one  specimen  the  average  mass  loss  for  each  test  Is 
considerably  higher  than  the  base  line  single  particle  value.  It  has  been 
postulated  that  this  Increase  In  average  mass  loss  Is  due  to  subsurface 
Interaction  of  damaged  material.  This  Increase  Is  frequently  referred  to 
as  "damage  enhancement",  and  a damage  enhancement  factor  has  been  defined 
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as  the  ratio  of  the  mass  loss  obtained  for  a pre*-damaged  surface  to  that 
obtained  for  a virgin  surface.  Damage  enhancement  Is  possibly  the  most 
sensitive  parameter  In  erosive  mass  loss.  Experimental  values  of  up  to 
3.3  have  been  recorded. It  Is  also  a complex  function  of  extrinsic 
variables  (such  as  velocity,  particle  size,  particle  type,  and  temperature), 
as  well  as  intrinsic  variables  (such  as  material  properties  and  material 
failure  modes). 


The  present  work  Is  directed  toward  quantifying  the  extent  and  nature 
of  the  subsurface  damage  produced  by  a variety  of  Impact  conditions  and  to 
assess  the  contribution  of  the  most  significant  aspects  of  the  subsurface 
damage  to  damage  enhancement.  The  ultimate  objective  will  be  to  suggest 
viable  micros true tural  modifications  that  would  Improve  the  erosion  re- 
sistance of  carbon  composites.  For  this  purpose,  detailed  post-test  exam- 
ination of  crater  formation  and  subsurface  damage  have  been  performed  for 
three  carbon-carbon  materials  (Table  1)  Impacted  by  glass,  nylon  or  water 
projectiles  In  the  velocity  range  1800  to  6000  ms  ^ at  temperature  between 
20  and  3600**C.  These  observations  have  been  correlated  with  Information 
obtained  from  high  speed  photography,  concerning  the  material  removal 
sequence,  and  from  selected  quasl-statlc  experiments, which  provide  ancillary 
data  concerning  both  the  conditions  for  damage  formation  and  the  physical 
properties  of  the  composite  constituents.  These  correlations  permit  per- 
llmlnary  postulates  to  be  constructed  about  the  details  of  the  cratering 
process  from  Initial  contact  to  damage  completion. 
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The  conditions  used  In  this  work  for  damage  analysis  are  primarily 
single  particle  Impacts  produced  by  1 mm  glass  beads  using  a high  energy 
capacitor  discharge  exploding  foil  system. The  test  facility  can 
perform  tests  at  temperatures  up  to  3600®C.  A commercial  Induction  heating 
unit  Is  used  to  heat  a 12  mm  diameter  by  12  mm  long  specimen  to  3600^0  In 
two  to  three  seconds. 

The  craters  produced  at  h3rpersonlc  particle  Impact  velocities  are 
examined  using  both  optical  and  electron  microscopy.  Initially  plan  views 
of  the  craters  are  obtained  using  scanning  electron  microscopy  to  provide 
a general  comparison  of  the  nature  of  the  damage  from  one  test  condition 
to  another  and  for  the  various  carbon-carbon  systems.  Then  the  crater  Is 
scanned  for  details  of  the  surface  fracture  modes.  Subsequently  the 
specimens  are  sectioned.  The  preferred  sectioning  method  Involves  Impreg- 
nation of  the  specimen  in  epoxy  to  preserve  the  details  of  the  microstructure 
during  subsequent  sectioning  operations  (the  epoxy  resin  can  be  made  to  com- 
pletely permeate  the  Interior  of  the  specimen  by  using  a prior  outgasslng 
procedure).  The  section  Is  prepared  using  a wire  saw  at  the  lowest  possible 
cutting  rate  and  then  metallographically  polished.  Optical  microscopic 
observations  are  made  on  progressive  cross  sections  to  establish  a three- 
dimensional  picture  of  the  extent  and  nature  of  the  subsurface  damage  and 
Its  relation  to  the  damage  observed  on  the  crater  surface.  The  general 
nature  of  the  damage  produced  for  a range  of  conditions  Is  summarized. 

2.1  DAMAGE  PERSPECTIVES 

A sequence  demonstrating  the  development  of  the  impact  damage  at  low 
temperatures  (0  - 2500“C)  as  a function  of  projectile  velocity  (for  1 mm 
glass  projectiles) Is  presented  in  Figure  1.  At  velocities  below  ~2000  ms 
the  projectile  is  retained  by  the  target  (Figure  la)  and  the  crater  consists 
primarily  of  a central  zone  that  accommodates  the  penetration  of  the  pro- 


6 


Rockwell  International 

Sctonce  Center 


SC5076.4TR 

jectile.  However,  at  higher  velocities  (Figure  lb),  the  projectile  is 
ejected,  and  the  crater  contains  an  outer  ledge  in  addition  to  a more 
deeply  penetrating  central  zone.  Surface  spallation,  manifested  as  the 
stripping  of  the  lateral  bundles  exposed  to  the  surface,  is  also  apparent 
(Figure  Ic)  ( longitudinal  bundle  Is  defined  as  a bundle  oriented  In  the 
direction  of  the  applied  force  , dynamic  or  quasi-static,  — usually  the 
bundles  normal  to  the  surface  — while  the  bundles  orthogonal  to  the  applied 
force  are  referred  to  as  lateral  bundles).  The  sample  depicted  In  Figure  Ic 
was  photographed  prior  to  Impact  to  demonstrate  that  the  damage  is  localized 
and  that  no  fine-scale  damage  occurs  beyond  the  area  encompassed  by  the 
fiber  stripping  boundary.  An  initial  evaluation  of  the  effect  of  micro- 
structure (see  Table  1)  on  the  damage  morphology  (Figure  Ic,  d,  e)  indicates 
that  as  the  proportion  of  longitudinal  bundles  increases,  the  central  pene- 
tration zone  tends  to  slightly  increase  in  depth,  while  the  outer  ledge 
diminishes  in  width.  Close  inspection  of  the  central  zone  (Figure  2a) 
indicates  the  existence  of  accommodation  damage,  manifested  as  exposed 
fibers  conforming  to  a spherical  boundary.  This  suggests  that  the  central 
zone  is  formed  by  the  direct  penetration  of  the  projectile  as  at  the  lower 
velocities.  Kinking  of  the  lateral  fiber  bundles  near  the  crater  periphery 
is  also  apparent  in  certain  microstructures  (Figure  2b),  and  the  presence 
of  short  lengths  of  longitudinal  fibers  laterally  displaced  away  from  the 
impact  center  (Figure  2c)  is  commonly  observed  on  longitudinal  bundles  at 
the  crater  ledge. 


The  damage  created  by  deformable  projectiles  (water  or  nylon)  has 
most  of  the  same  features  (Figure  3) , although  there  is  little  evidence  of 
accommodation  damage  in  the  central  zone.  Also,  under  multiple  impact 
conditions,  linking  of  the  outer  ledge  damage  zones  is  observed  at  closely 
spaced  Impacts. 

At  high  temperatures  (>  3000‘’C)  the  impact  damage  is  again  quite 
similar  (Figure  4),  except  that  the  impact  is  accompanied  by  permanent 
distortion  of  the  fibers  adjacent  to  the  central  zone  (Figure  4a)  and  by 
the  creation  of  permanent  interbundle  separations  (Figure  4b). 
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The  detailed  nature  of  the  impact  damage  has  been  obtained  from 
sectioned  specimens.  The  principal  modes  of  damage  identified  in  this 
study  are  summarized  in  schematic  overviews  (Figure  3a,  b,  c) . Sections 
close  to  the  center  of  impact  obtained  from  low  temperature  Impacts 
(<  3000®C)  show  that  the  longitudinal  bundles  directly  beneath  the  central 
penetration  zone  (Figure  5a)  are  extensively  kinked  (Figure  6)  and  that 
the  out-of-plane  lateral  bundles  are  vertically  displaced.  No  clear  re- 
lation between  the  kink  mode  (in  or  out-of-plane) , the  initiation  site,  or 
the  location  of  the  adjacent  out-of-plane  lateral  bundles  could  be  dis- 
cerned (Figure  6a).  Microfracture  is  invariably  observed  at  the  kink 
boundaries  (Figure  6a),  while  independent  macrofracture  (Figure  6b)  and 
partial  separation  of  the  kink  boundaries  (Flgu"  e 6c)  is  apparent  in  some 
longitudinal  bundles  near  the  crater  boundary.  Several  layers  of  in-plane 
lateral  bundles  (Figure  5b)  also  exhibit  kinking  that  extends  well  beyond 
the  periphery  of  the  crater  (Figure  7a) . These  kinks  are  mostly  in-plane 
(Figure  7b,  c)  and  generally  narrower  than  the  kinks  in  the  longitudinal 
bundles.  There  is  often  a distinct  fiber  orientation  relation  (Figure  7c), 
position  near  the  crater,  mixed  at  a distance  '7d^  from  the  center  of  im- 
pact (where  d^  is  the  bundle  thickness) , and  then  predominatly  negative  at 
larger  distances.  This  kinking  often  initiates  at  sites  close  to  the  comer 
of  out-of-plane  lateral  fiber  bundles.  A more  complex  mode  of  kinking  is 
apparent  in  the  lateral  bundles  directly  beneath  the  Impact  center  (Figure  8) . 
Inspection  of  these  fiber  failures  indicates  mixed  mode  kinking,  usually  on 
a plane  between  the  opposite  comers  of  the  adjacent  out-of-plane  lateral 
bundles  (Figure  8a).  However,  some  tensile  fractures  (Figure  8b)  and 
simple  kinking  adjacent  to  the  matrix  (Figure  8c)  can  also  be  detected. 

The  effects  of  varying  the  dimensions  of  the  longitudinal  bmdles  or 
the  projectile  velocity  on  the  extent  and  morphology  of  the  kinking  are 
summarized  in  Figures  9 and  10.  There  is  an  Increase  in  the  damage  extent 
with  the  projectile  velocity.  This  is  particularly  prominent  for  the  damage 
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In  the  lateral  fibers  (Figure  9) . However,  no  systematic  trend  In  the  ex- 
tent of  the  damage  zone  Is  apparent  for  the  various  bundle  configurations 
(Figure  9).  This  contrasts  with  some  significant  effects  of  the  bundle 
configuration  on  the  crater  dimensions  (Figure  1) . The  detailed  morphology 
of  the  kinking  In  tbe  longitudinal  bundles  does  vary  (Figure  10} . For 
relatively  large  bundle  diameters,  the  kinking  Is  confined  primarily  to  an 
Intensely  damaged  single  bundle,  with  some  partial  kinking  of  the  adjacent 
bundles.  For  smaller  bundles  the  kinking  Is  usually  less  Intense  and  dis- 
tributed over  two  or  more  adjacent  bundles. 


Near  the  crater  periphery  (Figure  5c)  lateral  bundle  fracture  Is 
typical  (Figure  11a,  b).  Occasionally,  sub?3urface  kinks  can  be  detected 
In  longitudinal  bundles  (Figure  11c)  at  the  lateral  bundle/matrix  Interface. 
It  Is  also  apparent  from  the  fiber  configuration  at  the  crater  surface 
(Figure  lid)  that  the  fractures  In  the  longitudinal  bundles  that  lead  to 
material  ejection  are  relatively  planar  In  character  and  usually  occur  In 
the  plane  of  the  lateral  bundle/matrix  Interface.  Some  Isolated  fractures 
In  the  longitudinal  bundles  are  also  evident. 


The  Impact  damage  at  elevated  temperatures  (>3000°C)  exhibit  some  Im- 
portant differences  from  those  discerned  at  low  temperatures.  Extensive 
permanent  distortion  of  the  bundles  near  the  crater  surface  (Figure  12) 
Indicate  that  the  fibers  are  capable  of  limited  plastic  extension  prior 
to  fracture.  The  kinking  Is  less  extensive,  especially  along  the  lateral 
bundles,  and  the  kinks  usually  arrest  before  they  propagate  through  the 
fiber  bundles  (Figure  12b).  Fiber  fracture  does  not  Invariably  occur  at 
the  kink  boundary  (Figure  12b) ; rather,  the  kink  strain  Is  accommodated. 

In  a significant  proportion  of  fibers,  by  a permanent  distortion  mechanism. 

2.3  MATERIAL  EJECTION 


High  speed  photographic  sequences  have  been  taken  of  a limited  number 
of  Impacts  on  2-2-3.  A sequence  obtained  for  a 1 mm  glass  projectile 
Impacting  at  2440  ms”^,  shown  In  Figure  13,  Indicates  the  existence  of  an 
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Initial  conical  spray  of  material  ejection,  followed  by  the  primary  ejection 
phase,  consisting  of  material  ejected  normal  to  the  specimen  surface  as  a 
central  plume  rising  out  of  the  crater.  (Note  that  the  diameter  of  the  base 
of  the  pltime  Is  approximately  equal  to  the  particle  diameter,  1 mm.)  By 
tracing  the  temporal  changes  In  the  location  of  a distinctive  feature  In 
the  plume,  the  approximate  time  at  which  It  was  freed  from  the  crater  can 
be  reconstructed  as  plotted  In  Figure  14.  It  Is  evident  that  although  the 
penetration  phase  was  completed  with  <20)is,  the  material  ejection  time  ex- 
ceeds 200ys — over  an  order  of  magnitude  longer.  The  primary  material 
ejection  thus  occurs  well  after  full  penetration. 
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Beside  direct  observations  of  the  Impact  process  and  post-test  exami- 
nations, quantitative  characterization  procedures  were  explored  In  order  to 
relate  the  dominant  failure  modes  to  material  properties.  In  general,  con- 
ventional test  procedures  are  not  directly  applicable;  the  strength  of 
the  composite  must  be  considered  on  a scale  relevant  to  the  size  of  the 
Impacting  particle.  The  tests  described  here  are  a step  In  that  direction. 
The  microscopic  failure  modes  associated  with  these  tests  are  compared  with 
those  found  for  particle  Impacts  as  described  In  Section  2.0. 


3.1  PENETRATION  TESTS 


An  Instrumented  drop  tower  was  used  for  the  penetration  tests^^^^. 
This  test  could  be  considered  dynamic  because  It  Involves  penetration 
velocities  of  a few  meters  per  second  (four  orders  of  magnitude  greater 
than  conventional  static  loading  rates) . These  tests  typically  last  a 
few  milliseconds,  which  Is  convenient  for  high  temperature  measurements. 

By  directly  Instrumenting  the  penetrator,  reasonably  accurate  simultaneous 
measurements  of  the  load  and  deflection  could  be  obtained.  A deflection 
gauge  system  was  used  for  penetration  measurements.  Two  Blomatlon  Model 
802  Transient  Records  were  used  to  record  load  and  deflection  signals.  The 
records  were  synchronized  by  a Tetronlx  FG  501  Function  Generator.  The 
Information  was  digitized,  fed  Into  a Hewlett-Packard  Model  9830A  computer, 
and  stored  on  magnetic  tape.  The  raw  data  could  then  be  plotted  or  further 
processed.  The  final  result  Is  a plot  of  penetration  pressure,  p,  versus 
depth  of  penetration  q/D  (normalized  to  the  Indenter  diameter) . A curve 
fitting  procedure  of  the  form 


P - (q/D)“  (1) 

Is  then  applied  to  the  data.  The  coefficient  p^  Is  the  pressure  at  one 
diameter  of  penetration,  while  the  exponent  n Is  the  slope  of  the  pressure/ 
penetration  curve.  The  parameters  p^  and  n are  determined  by  a least  mean 
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squares  fit.  An  example  of  this  analysis  Is  shown  In  Figure  15  for  a pene- 
tration test  on  2-2-3  material  at  IGOO^C.  Penetration  pressures,  p^, 
obtained  on  2-2-3  carbon-carbon  at  temperatures  up  to  3600‘’C,  are  plotted 
In  Figure  16;  the  temperature  dependence  of  the  mass  loss  Is  also  plotted 
for  comparison.  Comparing  the  form  of  the  mass  loss  data  with  the  penetra- 
tion test  results  Indicates  that  the  correlation  between  the  two  damage  modes 
Is  not  direct,  but  additional  thermomechanical  parameters  are  required  to 
establish  an  accurate  model  of  the  mass  removal.  These  relations  will  be  ex- 
plored further. 


The  damage  observed  In  the  penetration  samples  Is  mostly  similar  In 
character  to  that  obtained  on  the  Impact  samples,  but  the  scale  of  the 
damage  Is  larger  (because  of  the  larger  contact  area) . The  extended  scale 
has  the  Important  advantage  that  the  damage  sequences  can  be  more  effectively 
traced,  and  this  provides  some  Insights  into  the  origins  of  the  Impact  damage. 
The  damage  In  the  longitudinal  bundles  has  several  Important  features  (Figure 
17a)  that  confirm  Inferences  made  from  the  damage  In  the  particle  Impact 
samples.  The  kinking  In  the  bundles  beneath  the  Impact  crater  (bundles  c 
and  d,  and  b,  e,  f below  lateral  bundle  3)  appears  to  exhibit  little  spatial 
correlation,  either  along  the  bundles  or  between  adjacent  bundles.  Rink  pairs 
of  opposite  sign  appear  to  be  the  primary  mechanism  for  retaining  the  Initial 
bundle  alignment.  Conversely,  the  kinking  In  the  bundles  beside  the  crater 
(bundles  b and  f)  tend  to  exhibit  a uniform  spatial  separation, with  the 
kinks  usually  Initiating  near  the  comer  of  an  out-of-plane  lateral  bundle 
(e.g.,  lb,  2b,  If,  2f,  3f ) . The  kink  separation  Is  thus  approximately 
equal  to  the  diameter  of  a lateral  bundle  plus  a matrix  Interstice.  Also, 
the  kinks  tend  to  conform  to  a single  orientation.  Finally,  It  Is  noted 
that  kinks  occasionally  occur  In  bundles  that  are  not  directly  subject  to 
Interaction  with  the  penetration  (e.g.,  bundle  a). 


The  damage  to  the  In-plane  lateral  bundles  also  has  Important  impli- 
cations (Figure  17b).  The  kinks  In  the  bundles  beside  the  crater  are  quite 
restricted,  in  marked  constrast  to  the  extensive  kinking  In  the  high  velocity 
impact  samples.  These  kinks  are  highly  correlated.  Initiating  from  the 
comers  of  out-of-plane  lateral  fiber  bundles  along  a single  plane.  It  is 
Instmctlve  to  regard  these  kinks  as  being  Induced  by  a shear  displacement, 
as  depicted  by  the  vectors  in  Figure  17b. 
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The  compliance  of  each  zone  of  the  composite  has  been  evaluated  using 
a mlcTOCompresslve  test.^^^^  The  force  deflection  behavior  Indicates  that 
the  longitudinal  bundles  exhibit  the  maximum  stiffness  E (=  30  GPa) , i^lle 


the  lateral  fibers  (In  both  locations)  and  the  matrix  have  a similar 

stiffness  (E  3 5 GPa)  at  low  stress  levels.  The  greater  stiffness  of  fiber 

composites  In  the  longitudinal  direction  Is  a well-documented  effect;  the 

ratio  of  the  lateral  to  longitudinal  stiffness  depends  on  the  fiber  loading, 

(12) 

the  matrix  material,  etc.  However,  the  matrix  becomes  nonlinear  at 

larger  stresses  (due  to  a permanent  compression  associated  with  the  porosity) , 
and  thereafter  the  effective  tangent  modulus  decreases  continuously  with 
strain.  The  deformation  characteristics  of  the  matrix  beyond  yield  can  be 
assessed  from  a spherical  Indentation  test,  using  an  analysis  comparable 
to  that  described  for  the  penetration  tests  (Section  3.1).  Then,  by  noting 
that  the  Indentation  pressure  p Is  '3o^  (where  o^  Is  the  uniaxial  flow 
stress),  the  test  Indicates  that  the  Initial  flow  stress  of  the  matrix  Is 
-30  MPa. 

3.3  PLANE-STRAIN  COMPRESSION  TESTS 

Quantitative  evaluation  of  the  stresses  required  for  kink  formation 
has  been  achieved  using  a plane-strain  compression  test  of  the  type  Illus- 
trated schematically  In  Figure  18.  The  test  entails  compressing  a rec- 
tangular sample  (approximate  dimensions  25  x 5 x 5 mm)  between  two  parallel 
plattens  (5  mm  In  width)  and  applying  orthogonal  constraints  to  suppress 
fracture.  Parallelism  of  the  plattens  was  achieved  using  a conventional 
hemispherical  alignment  fixture.  The  constraint  was  applied  using  horizontal 
mlnlclamps  and  the  forces  measured  using  strain  gauges.  Compression  was 
then  applied,  at  a constant  rate  of  displacement,  using  an  Instron  testing 
machine.  A typical  load  displacement  curve  is  plotted  In  Figure  18,  in- 
dicating a load  Instability  prior  to  the  onset  of  nonlinearity. 
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Saaples  tested  to  pre-selected  loads  were  subsequently  sectioned  and 
metallographlcally  polished.  No  permanent  damage  could  be  detected  In  the 
linear  region,  but  kinks  In  one  or  more  of  the  longitudinal  fibers  Invariably 
could  be  detected  after  the  load  Instability  (Figure  18).  The  Instability 
Is  thus  presumed  to  coincide  with  the  onset  of  kinking  In  the  longitudinal 
fiber  bundles.  The  kinks  exhibit  the  same  morphology  (Figure  19)  as  those 
observed  In  the  Impact  or  penetration  tests.  At  larger  loads  (Figure  18), 
compressions!  failures  of  the  lateral  bundles  could  be  detected  In  the  form 
of  complex  kinks  between  out-of-plane  lateral  bundles,  as  observed  In  the 
Impacted  samples  (Figure  8). 


The  stresses  In  the  bundles  at  the  damage  thresholds  were  computed 
by  taking  the  applied  threshold  displacements  and  converting  them  Into 
strains  and  then  stresses,  using  the  moduli  (for  the  longitudinal  or 
lateral  bundles,  or  the  matrix)  computed  In  the  compliance  tests  (Section 
3.2).  The  validity  of  the  computed  stress  levels  was  approximately  sub- 
stantiated by  comparing  the  threshold  stresses  In  the  longitudinal  fibers 
with  those  obtained  by  assuming  that  the  load  up  to  Instability  Is  trans- 
mitted exclusively  by  the  longitudinal  bundles  (because  of  their  much 
larger  modulus).  The  threshold  stresses  derived  In  this  manner  for 
2-2-3  are  plotted  In  Figure  20  as  a function  of  the  constraining  stresses. 
The  threshold  stress  for  longitudinal  kinking  Is  larger,  even  though  It 
occurs  at  a lower  applied  force  since  the  longitudinal  bundles  experience 
a substantially  larger  stress. 


The  Increase  In  the  threshold  stress  for  longitudinal  bundle  kinking 

as  the  constraint  Is  Increased  Is  approximately  proportional  to  the  level 

of  the  constraint.  The  result  suggests  a kink  Initiation  criterion  dictated 

by  the  maximum  shear  stress  (l.e.,  the  difference  between  the  maximum  and 

minimum  principal  stresses),  as  also  Inferred  In  a recent  detailed  study 

of  the  kink  formation  process. In  addition,  under  dynamic  conditions 

there  might  also  be  a minimum  stress  duration  for  kink  Initiation:  com- 

0 3) 

parable  to  the  Incubation  time  for  crack  activation. 
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The  laltlatlon  of  compresslonal  failure  In  the  lateral  bundles  does 


not  exhibit  a strong  dependence  on  the  applied  constraint,  suggesting  a 
critical  normal  stress  criterion  as  an  approximate  descriptor  of  damage 


Initiation.  This  Is  again  consistent  with  an  Independent  analysis  of 

(14) 

the  compresslonal  damage  process. 
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The  extent  of  the  damage  to  carbon  composites  created  by  projectile 
lsq>act  Is  dictated  by  the  stresses  that  develop  due  to  the  Impact  event 
and  the  material's  resistance  to  the  important  modes  of  damage.  Since 
damage  can  derive  from  both  the  compressive  stresses  (kinking)  and  tensile 
stresses  (fracture),  some  Idea  of  the  general  nature  of  the  stresses  In- 
volved In  the  entire  Impact  cycle  Is  needed  to  Interpret  the  damage  sequence. 
Then,  knowledge  of  the  conditions  for  Initiating  the  various  modes  of  damage 
Is  required  In  order  to  correlate  the  dynamic  stresses  with  both  the  In- 
cidence of  damage  and  material  removal.  A detailed  analysis  of  this  complex 
process  can  only  be  performed  using  numerical  methods.  However,  a per- 
spective of  the  process  can  be  obtained  from  simplified  analyses  which 
allow  elucidation  of  several  Important  principles  for  material  Improvement 
with  respect  to  erosive  environments. 

4.1  IMPACT  STRESSES 


FoUotrlng  the  Initial  contact  of  the  projectile  and  the  target,  a 
shock  wave  Initiates  and  propagates  away  from  the  contact  zone  In  both 
the  target  and  the  projectile.  Due  to  the  localized  nature  of  the  Impact, 
the  essentially  spherical  shock  pulse  attenuates  and  waves  are  reflected 
at  the  existing  free  surface  and  those  produced  during  the  penetration  pro- 
cess. The  cummulatlve  wave  Interactions  In  the  projectile  produce  significant 
tensile  stresses.  These  eventually  develop  along  the  normal  to  the  Inter- 
face and  the  mutual  contact  terminates.  The  tensile  stresses  In  the 
projectile  also  cause  the  fracture  and  fragmentation  observed  In  brittle 
projectiles  such  as  glass.  The  target  material  at  the  shock  front  Is 
subjected  to  a hydrostatic  compression. 


The  shock  reflections  at  the  target  surface  produce  near-surface 
rarefaction  waves  (Figure  21) . These  reflected  waves  are  the  probable 
cause  of  the  Initial  spallation  damage  (Figure  Ic)  occurring  on  a time 
scale  dictated  by  the  velocity  of  the  attenuating  shock  wave.  The 
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roughened  surface  of  Che  contact  zone  produced  by  spallation  will  sub- 
sequently be  enveloped  by  the  sphere  and  carried  downward  during  Che 
sphere  penetration  stage.  The  presence  of  this  roughened  zone  on  the 
contact  Interface  could  produce  minor  perturbations  In  Che  early  target/ 
sphere  Interfaclal  pressure  distribution.  The  sphere  continues  to  deform 
as  It  penetrates,  so  Che  particle  velocity  and  pressure  at  the  actual 
carget/sphere  Interface  are  fairly  complex  functions  of  the  depth  of  pene- 
tration. Also  Che  rarefaction  wave  developed  by  reflection  from  the 
target  surface  (as  the  shock  wave  moves  outward)  continues  to  develop^ but 
Is  complicated  by  Che  presence  of  Che  crater  which  Introduces  a second 

I free  surface. 

I 

I 

I 


The  shock  related  phenomena  comprise  only  Che  very  earliest  stages 
of  the  penetration  process.  The  expanding  shock  front  quickly  attenuates 
and  rapidly  moves  out  of  the  primary  damage  zone.  As  penetration  continues 
the  shock  wave  Is  succeeded  by  elastic/plastic  waves  (Figure  21) . The 
term  plastic  wave  is  used  here  to  denote  the  relatively  slowly  moving 
waves  Chat  encompass  Che  non-linear  zone  of  material  response.  The  non- 
llnearlcles  could  be  attributed  to  matrix  densification,  matrix  or  fiber 
plasticity,  or  viscoelastic  effects.  The  elastic  wave  generates  radial 
compression  and  tangential  tension  with  the  approximate  relative  magnitudes 
depleted  In  Figure  21;  while  Che  plastic  wave  only  generates  compressive 
stresses. with  the  radial  compression  substantially  larger  than  the  tan- 
gential compression  (the  difference  being  dictated  by  the  dynamic  plastic 
flow  stress). 


The  penetration  proceeds  as  Che  projectile  decelerates,  even  Chough 
it  may  be  highly  fractured  or  extensively  deformed.  This  process,  coupled 
with  the  Inertia  of  the  expanding  boundary,  continues  to  impose  a direct 
pressure  on  the  surrounding  material.  The  plastic  zone  continuously  in- 
creases In  extent  during  this  phase. 
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When  penetration  terminates,  the  elastic  unloading  waves  become 
fully  developed.  These  waves  can  Interact  with  the  plastic  waves  to  generate 
substantial  tensile  stresses  In  orientations  approximately  normal  to  the 
plastic  wave  front. Material  removal  Is  typically  attributed  to  the 
large  tensile  stresses  that  develop  after  full  penetration. 


These  effects  have  been  described  for  a target  with  Isotropic  physical 
properties.  Local  perturbations  are  likely  for  carbon-carbon  composites. 
These  would  substantially  modify  the  details  of  the  stress  magnitude  and 
orientation  (e.g.,  orienting  the  normal  stresses  along  the  bundle  axis), 
but  the  sequences  and  modes  of  stress  development  should  be  essentially 
unchanged. 


4.2  DAMAGE  MECHANISMS 

The  character  of  the  residual  damage  In  Impacted  samples,  combined 
with  the  observed  material  ejection  sequence  and  the  likely  states  of 
stress  developed  during  Impact,  suggest  the  following  evolution  of  damage. 

The  shock  wave  that  develops  after  Initial  contact  probably  leads  to  the 
formation  of  surface  spallation  (fiber  stripping)  due  to  rarefaction 
waves  that  derive  from  free  surface  reflections.  The  large  compressive 
stresses  at  the  shock  front  might  also  generate  compresslonal  failure  of 
both  the  lateral  bundles  beneath  the  Impact  center  and  the  lateral  or  longi- 
tudinal bundles  adjacent  to  the  surface  (this  process  could  also  contribute 
to  the  attenuation  of  the  shock  wave).  However,  kinking  may  not 
be  associated  with  the  Initial  passage  of  the  shock  front.  If  the 
shear  stresses  are  small.  During  subsequent  penetration,  when  the  shock 
condition  Is  succeeded  by  elastlc/plastlc  waves,  compresslonal  and  shear 
components  of  stress  can  exist  and  the  damage  consists  of  both  compresslonal 
bundle  failure  and  kinking.  The  extentlon  of  these  damage  zones  is  dictated 
by  the  spatial  and  temporal  details  of  the  dynamic  stress  field.  The  damage 
processes  are  Illustrated  schematically  (Figure  22) , using  possible  variations 
in  the  dynamic  stresses  along  the  penetration  (z)  axis,  with  critical  damage 
stresses  superimposed.  The  compresslonal  failure  zone  and  the  kinking 
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zone  z^  , are  located  by  recalling  that  the  stress  (and  wave  speed)  trans- 
max 

mltted  to  the  lateral  bundles  Is  lover  than  that  transmitted  along  the 
longitudinal  bundles  (because  of  differential  modulus  effects)  and  that 
the  stress  must  be  applied  for  a minimum  (incubation)  time. 


When  penetration  is  complete,  the  tensile  stresses  that  develop  due 
to  interaction  of  the  unloading  waves  with  the  plastic  waves  are  presumed 
to  cause  a portion  of  the  damage  zone  to  be  ejected.  The  extent  of  the 
ejected  zone  depends  on  the  residual  tensile  strength  of  the  kinked 
material  (In  orientations  normal  to  the  tensile  stress):  a property  that 
has  not  yet  been  adequately  studied  (see  Figure  6c).  This  postulate  is 
consistent  with  the  primary  material  removal  occurring  well  after  full 
penetration  and  with  the  ejecta  trajectories  being  approximately  normal 
to  the  target  surface  (Figure  13) . 


5.0  SUMMARY  AND  GENERAL  IMPLICATIONS 
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The  considerable  differences  in  the  extent  and  character  of  the 
kinking  damage  beneath  and  beside  the  penetration  zone  have  Important 
implications  for  the  Interpretation  of  the  damage  process  and  for  de- 
signing more  damage  resistant  materials. 


Kinking  of  the  longitudinal  fiber  bundles  beneath  the  penetration  zone 
is  reproduced  under  djmamic  and  quasi-static  conditions;  its  extent  is 
probably  dictated  by  the  stress  state  near  the  final  stages  of  penetration 
when  the  projectile  is  moving  relatively  slowly  and  the  applied  pressure 
is  comparable  to  the  quasi-static  penetration  pressure.  The  effects  of 
the  bundle  dimensions  and  the  matrix  and  fiber  properties  on  the  resulting 
impact  damage  should  thus  be  related  to  their  separate  effects  on  the 
critical  stress  for  kinking  and  the  penetration  pressure  — the  latter 
through  its  influence  on  both  the  penetration  distance  and  the  dynamic 
stresses.  It  is  already  apparent  that  kinking  is  suppressed  by  increasing 
the  matrix  yield  strength  and  elastic  modulus (an  immediate  Implication 
for  materials  development) . Whereas  the  observed  lack  of  correlation  of 
kinks  in  adjacent  bundles  suggests  that  the  influence  of  the  bundle  di- 
mensions on  the  kink  formation  condition  is  small.  Similarly,  the  pene- 
tration pressure  should  be  enhanced  by  increasing  the  matrix  flow  strength; 
while  the  effect  of  the  bundle  dimensions,  although  uncertain,  is  probably 
quite  small.  However,  it  is  not  immediately  obvious  that  the  extent  of 
the  kinked  zone  is  decreased  by  increasing  the  penetration  pressure.  This 
ambiguity  arises  because  the  diminution  in  penetration  distance  is 
counteracted  by  an  increase  in  the  amplitude  of  the  d5mamic  stresses.  Further 
work  is  clearly  needed  to  develop  a detailed  picture  of  the  role  of  the 
penetration  pressure  on  the  extent  of  longitudinal  fiber  bundle  kinking  and 
its  dependence  on  microstructure. 


The  kinking  of  the  lateral  fiber  bundles  beside  the  penetration  zone  is 
only  detected  under  fully  dynamic  conditions,  and  its  extent  is  strongly 
dependent  on  the  impact  velocity.  Presumably,  therefore,  the  conditions 
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that  encourage  this  kinking  develop  earlier  In  the  penetration  process; 
perhaps  because  the  pressure  applied  to  the  bundles  declines  before  pene- 
tration Is  complete.  The  morphology  of  the  kinks  also  suggests  that  the 
free  surface  has  an  Important  Influence  on  the  kinking  process. 

The  longitudinal  fiber  bundles  beside  the  penetration  zone  must  also 
be  damaged  during  the  Impact  cycle  to  account  for  the  observed  occurrence 
of  the  outer  ledge  (Figure  lb).  Since  there  is  no  systematic  evidence  of 
residual  damage  in  these  bundles,  suggestions  concerning  the  nature  of  this 
damage  can  only  be  speculative.  There  is  much  supportive  evidence  for  the 
speculation  that  these  bundles  are  damaged  by  kinking  during  the  penetration 
phase  and  that  the  severity  of  kinking  leads  to  the  subsequent,  full  ejection 
of  the  damaged  zone.  The  regularly  spaced  kinks  (that  initiate  at  the  matrix/ 
bundle  interface)  observed  adjacent  to  the  impact  center  in  the  penetration 
tests,  indicate  that  the  relative  displacements  that  accompany  penetration 
can  generate  in-plane  shear  stresses,  as  indicated  schematically  in  Figure 
23.  Then,  the  sequence  of  kinks  required  to  accommodate  an  applied  in-plane 
shear  stress,  when  the  restraint  provided  by  the  orthogonal  bundles  has  been 
removed  by  prior  kinking,  suggests  that  a series  of  shear  kinks  could  de- 
velop beside  the  penetration  zone.  In  support  of  this  hypothesis  we  note 
that  the  relatively  planar  character  of  the  fractured  longitudinal  bundles 
at  the  base  of  the  outer  ledge  is  typical  of  kink  boundary  fracture,  while 
the  common  occurrence  of  these  bundle  fractures  at  a lateral  bundle/matrix 
interface  is  typical  of  shear  induced  kinks.  Also,  the  short  fiber  segments 
sometimes  observed  on  the  longitudinal  bundles.  Figure  2c,  could  be  the 
remnants  of  the  kinked  zone.  Should  this  damage  process  pertain,  the  extent 
of  the  damage  normal  to  the  impact  direction  would  exhibit  (in  addition  to 
the  prevading  influence  of  the  matrix  flow  stress  and  modulus)  a significant 
dependence  on  the  proportion  of  longitudinal  fibers^ in  the  sense  that  the 
damage  zone  will  decrease  as  the  content  of  longitudinal  fibers  increases. 

This  effect  is  consistent  with  the  observations  of  the  effect  of  the  longi- 
tudinal fiber  content  on  the  extent  of  the  outer  ledge.  It  is  also  noted, 
however,  that  kinking  of  the  restraining  lateral  bundles  is  a necessary  pre- 
cursor to  the  damage  process;  this  requirement  implies  an  upper  limit  on 
the  zone  size  that  might  be  encountered  for  relatively  small  longitudinal 
fiber  contents. 
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Figure  1 


Figure  2 


Figure  3 
Figure  4 


Figure  5 


Figure  6 


Figure  7 


Scanning  electron  microscope  overview  of  craters  formed  by 
1 mm  glass  projectile  impacts. 

Scanning  electron  micrographs  of  details  of  impact  craters 
in  2-2-3  material  (a)  accommodation  damage  at  the  base  of  the 
central  zone  (b)  kinking  of  a lateral  bundle  at  the  crater 
edge  (c)  short  fiber  segments  on  a longitudinal  bundle  at  the 
outer  ledge. 

Impact  damage  on  2-2-3  created  by  1 mm  v/ater  drops  at  4000  ms~^. 
Scanning  electron  micrographs  of  high  temperature  (3500®C) 
impact  damage  in  1-1-5  material  impacted  by  a 1 mm  glass 
projectile  at  6000  ms~^  indicating  (a)  distortion  of  longitudinal 
fibers  in  the  central  zone,  (b)  separation  betv/een  the  lateral 
bundles  and  the  matrix  at  the  outer  edge. 

Schematics  of  the  damage  created  by  projectile  impact  (a)  a 
section  near  the  impact  center  revealing  longitudinal  bundles, 

(b)  a section  near  the  impact  center  revealing  lateral  bundles, 

(c)  a section  near  the  crater  periphery. 

Optical  micrographs  of  polished  sections  near  the  impact  center 
revealing  longitudinal  bundles  (a)  kinking  damage  beneath  the 
center  of  impact,  (b)  bundle  fracture  beside  the  crater,  (c) 
partial  separation  of  bundle  at  kink  boundaries. 

Optical  micrographs  of  polished  sections  near  the  impact  center 
revealing  lateral  bundles  (a)  an  overview  of  the  kink  locations, 

(b)  a detail  showing  the  local  consistency  of  the  kink  morphology, 

(c)  a detail  showing  a kink  beside  the  crater  wall. 
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Figure  8 


Figure  9 


Figure  10 


Figure  11 


Figure  12 
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Optical  micrographs  of  a section  near  the  impact  center  shov/ing 
compressional  damage  in  lateral  bundles  beneath  the  impact 
zone  (a)  mixed  mode  kinks  adjacent  to  out-of-plane  lateral 
bundles,  (b)  distortion  at  the  crater  base,  (c)  simple  kinking 
and  fracture  remote  from  the  crater. 

The  maximum  extent  of  the  kinking  in  the  longitudinal  and  lateral 
bundles  as  a function  of  the  projectile  (1  mm  glass  sphere) 
velocity  for  three  microstructures. 

Optical  micrographs  of  sections  near  the  impact  center  revealing 
longitudinal  bundles,  illustrating  the  effects  of  bundle 
dimensions  on  the  kink  morphology:  1 mm  glass  sphere  impacts 
at  6000  ms”\ 

Micrographs  of  fractures  near  the  crater  perphery  (a)  a 
scanning  electron  micrograph  of  a fractured  lateral  bundle, 

(b)  an  optical  micrograph  of  a polished  section  showing  a 
lateral  bundle  fracture,  (c)  a scanning  electron  micrograph 
of  a polished  section  showing  a sub-surface  kink  in  a 
longitudinal  bundle,  (d)  an  optical  micrograph  of  a polished 
section  showing  the  morphology  of  the  longitudinal  bundle  fractures. 
Optical  micrographs  of  sections  through  a 2-2-3  material 
impacted  at  3500*C  by  a glass  projectile  Jit  6000  ms"^  (a)  an 
overview  showing  lateral  bundle  distortion,  (b)  an  overview 
of  longitudinal  bundle  kinking,  (d)  detail  of  a kink  in  a 


longitudinal  bundle. 

Figure  13  A high  speed  photographic  sequence  of  a 1 mm  glass  sphere  impact 
at  2400  ms"^ . 
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Figure  14 

Figure  15 

Figure  16 

Figure  17 


Figure  18 

Figure  19 


Figure  20 

Figure  21 


Figure  22 


Figure  23 


The  temporal  development  of  the  processes  accompanying  a 
1 mm  glass  sphere  impact  at  2440  ms”\ 

The  penetration  pressure  for  2-2-3  material  at  lOOO^C,  obtained 
using  an  instrumented  drop  weight  test. 

The  temperature  dependencies  of  the  penetration  pressure,  p^, 
and  the  mass  loss  ratio  for  2-2-3  material. 

Optical  micrograph  of  a polished  section  from  a drop  weight 
sample  of  2-2-3  revealing  (a)  the  longitudinal  bundles, 

(b)  the  lateral  bundles  beside  the  impact  center. 

A schematic  of  the  plane  strain  compression  test,  and  a 
typical  load,  deflection  curve. 

An  optical  micrograph  of  a polished  section  through  a compressed 
plane  strain  test  sample  revealing  kinks  in  the  longitudinal 
bundles. 

The  threshold  stresses  for  damage  formation  plotted  as  a function 
of  the  lateral  constraint. 

A schematic  indicating  a possible  sequence  of -dynamic  stress 
development  during  hypervelocity  impact  in  graphite  based 
materials. 

A schematic  showing  the  possible  temporal  and  spatial  character 
of  the  principal  and  shear  stress  development  along  the  penetration 
axis,  during  the  projectile  penetration  phase,  with  critical 
damage  formation  stresses  superimposed. 

A schematic  indicating  the  possible  development  of  kinks  in 
longitudinal  bundles  at  the  outer  ledge. 
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One  of  the  important  modes  of  compressions!  damage  in  carbon 
composites  - complex  kinking  of  the  lateral  bundles  - has  been  studied. 

It  is  shown  that  damage  propagation  is  a consequence  of  matrix  plasticity, 
with  the  damage  trajectory  approximately  coinciding  with  the 
elastic/plastic  boundary.  The  analysis  indicates  that  the  damage 
resistance  is  enhanced  by  increasing  the  matrix  yield  strength  and  by 
optimizing  the  relative  dimensions  of  the  bundles. 
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Recent  studies  have  indicated  the  existence  of  at  least  two  modes  of 
compressional  damage  in  three-dimensional  composites. The  more 
extensive  damage  mode  is  bundle  kinking:  the  characteristics  of  which 
have  been  explored  in  a previous  study. 1 The  second  mode  is  a complex 
fracture  process^  that  occurs  due  to  the  compression  of  a lateral  bundle 
between  (two  opposing)  orthogonal  lateral  bundles.  This  damage  mode  is 
referred  to  as  complex  kinking  (for  reasons  which  will  become  apparent 
later) . 

The  objective  of  this  study  is  to  identify  the  primary  mechanism  of 
complex  kinking  and  thereafter,  to  deduce  the  microstructural  parameters 
that  affect  its  incidence.  The  study  is  part  of  a larger  program  devoted 
to  the  comprehensive  characterization  of  the  damage  modes  in  carbon 
composites,  as  required  to  design  microstructures  with  an  enhanced 
resistance  to  erosion  or  structural  degradation. ^ 

2.  DAMAGE  OBSERVATIONS 

Observations  of  lateral  fiber  damage  have  been  made  in  carbon 
composites  after  either  impact  by  small  projectiles  at  high  velocities,  or 
quasi-static  plane  strain  compression. ^ Some  typical  damage  zones  are 
depicted  in  Fig.  1.  The  damage  usually  (Fig.  la),  but  not  always 
(Fig.  lb),  occupies  a trajectory  between  opposite  corners  of  the 
orthogonal  lateral  bundles.  The  damage  zone  is  usually  narrow,  and 
contains  short  segments  of  fractured  fibers.  The  fractured  fibers 
frequently  exhibit  two  preferred  residual  orientations,  especially  when 
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the  zone  width  exceeds  the  fracture  length  ( A in  Fig  .la  ); 
otherwise,  there  are  no  systematic  orientation  effects.  The  orientation 
relations,  when  they  exist,  are  reminiscent  of  those  that  prevail  for 
kinking. 1 This  damage  mode  is  thus  referred  to  as  complex  kinking. 


3.  DAMAGE  ANALYSIS 

The  bases  for  quantitative  damage  analysis  arei,  a)  the  tensile 
stresses  that  develop  in  the  fibers,  and  b)  the  conditions  for  fiber 
fracture.  The  intent  of  the  present  analysis  is  to  develop  approximate 
stress  and  fracture  formulations,  using  simplifying  assumptions  where 
required,  that  permit  the  identification  of  the  damage  mechanism.  This 
analysis  allows  several  of  the  important  parameters  that  control  damage 
formation  to  be  elucidated,  and  forms  the  basis  for  a subsequent,  more 
rigorous,  damage  analysis. 

Initially,  the  forces  transmitted  to  the  bundles  must  be  related  to 
the  remote  compression,  These  forces  depend  on  the  properties  of  the 

composite  constituents.  Since  the  overall  compliance  of  the  composite  is 
dictated  primarily  by  the  constituent  with  the  largest  modulus,  i.e.  the 
longitudinal  bundles  (at  least  until  these  bundles  are  damaged  by 
kinking),  the  stress  transmitted  to  the  lateral  bundles  is  given  by: 

2 

o ” (^lat'^^long)  (^■'^lat'^'^long  ) 
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where  d is  the  bundle  diameter,  E is  Young's  modulus  and  the  subscripts 
lat.  and  long,  refer  to  the  lateral  and  longitudinal  bundles, 
respectively.  At  relatively  large  stress  levels  f beyond  the  interstice 
yield  stress } a tends  to  localize  at  the  orthogonal  lateral 

bundles.  However,  if  the  restraint  provided  by  the  longitudinal  bundles 
is  removed  by  kink  initiation,  the  stress  on  the  lateral  bundles  increases 
to. 


0 (1  + d,  ^/d-  )' 

“>  lat  long 


(2) 


The  response  of  the  fibers  in  the  lateral  bundles  to  the  localized 
stress  (7  dictates  the  damage.  Analysis  of  the  initial,  elastic  response 
is  extremely  complex  if  the  elastic  anistropies  of  the  bundles  are 
rigorously  included.  An  approximate  approach  is  thus  preferred  for 
present  purposes.  The  simplest  approximation  that  provides  a reasonable 
insight  into  stress  development  within  the  fibers  (Fig.  2)  involves  the 
estimation  of  the  displacements  in  the  bundle,  by  assuming  elastic 
, isotropy , and  then  uses  the  corresponding  fiber  radii  of  curvature  to 

' calculate  the  fiber  stresses,  by  neglecting  the  constraint  of  the  matrix. 

! 

The  regions  of  maximum  bundle  displacement  are  of  primary  concern. 

i These  occur  near  the  bundle  surface  (Fig.  2),  and  can  thus  be  obtained 

I with  sufficient  accuracy  by  employing  solutions  for  semi-infinite  bodies 

^ subjected  to  surface  tractions.  The  displacements  in  the  horizontal 

I direction  are  approximately  uniform  (Appendix)  and  do  not,  therefore, 

( 

I 

I 

i 

i 

I 
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contribute  to  fiber  curvature.  The  approximate  vertical  displacementi at 
the  surface,  VqCx),  are; 


"oW  “ He. 


. i J 


+ x Zn(z*/x)  + d. 


■C1-.W. 


(3) 


(x  < dj^/2) 


where  v is  Poisson's  ratio,  is  the  dimension  of  the  pressure  zone, 

X*  is  the  position  of  zero  displacement  (taken  in  this  instance  to  be 
half  the  in~plane  dimension  of  the  bundle),  and  <E>  is  an  average  modulus 
for  the  bundle.  The  radii  of  curvature  R(x)  of  the  fibers  are3- 


1 

RTx) 


d Vq(x) 

T1 

CiX 


2o  [l-2(x/dj^)] 
7r<E>x[l-  (x/dj^)] 


(4) 


The  tensile  stresses  in  the  fibers  are  thus  given  approximately  by; 


o Er  rr  / R 

max  i // 


2a(Ey./<E>)r^  Cl-2(x/d  )] 


TTx  [l-(x/d^)] 


(5) 


where  ia  the  maximum  tensile  stress  at  the  fiber  surface,  E^.  is  the 

fiber  modulus  and  rj  the  fiber  radius.  The  stress  is  thus  zero  of  the 
center  and  exhibits  a singularity  at  the  edge.  The  stress  near  the 
singularity  will  undoubtedly  be  modified  by  various  relaxation  effects; 
but  w can  reasonably  conclude  that  the  stress  in  the  fibers  near  the  edge 
of  the  loading  zone  will  be  large  enough  to  induce  fiber  fracture  at 
relatively  small  . The  extent  to  which  the  fracture  propagates  through 
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Che  bundle  is  dictated  by  the  stress  gradient.  Since  Che  stresses  are  a 
maximum  along  the  z-axis , the  displacements  v(x,z)  (Appendix)  at  small  x 
and  z are  of  interest.  The  second  derivative  of  these  displacements 
yields  Che  following  expression  for  Che  maximum  sub-surface  fiber  stress 
distribution; 


2o  \l^/<E>)Vf  p (6) 

°max  - (l+v)j 

The  principal  feature  to  note  is  that  the  stress  diminshes  rapidly  with 

distance  from  the  surface,  as  ~ reducing  to  values  of  the  order  of 

the  applied  pressure,  a , at  5-10  fiber  diameters.  It  seems  improbable 

00 

therefore  that  the  damage  propagation  could  extend  beyond  the  first  few 
fibers . 

We  will  thus  presume  that  matrix  plasticity  precedes  substantial  fiber 
fracture,  and  determine  whether  the  consequences  of  plasticity  are 
consistent  with  the  damage  observations.  The  zone  of  matrix  deformation 
that  accompanies  Che  onset  of  plasticity  can  be  approximately  assessed 
from  slip  line  field  solutions. 4 (it  should  be  noted,  however,  that 
this  mode  of  analysis  implicitly  assumes  a plastic/rigid  continuum  and 
thus  neglects  elastic  effects).  Slip  line  fields  derived  for  aligned  and 
misaligned  compression  between  parallel  platCens  are  plotted  in  Fig.  3. 
Marked  similarities  between  the  outer  slip  line  trajectories 
(AOb' ,A'0B)and  Che  damage  zones  (Fig.  la)  are  immediately  apparent. 

These  outer  trajectories  represent  the  plastic  boundary,  across  which  the 
tensile  component  of  velocity  is  discontinuous.  This  discontinuity 


i 


i 
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corresponds  to  a zero  radius  of  curvature  for  the  fibers  and  hence,  an 
infinite  tensile  stress.  The  magnitude  of  the  fiber  stress  at  the  plastic 
boundary  will,  of  course,  be  modified  when  elastic  effects  are 
incorporated;  but  the  stresses  are  still  expected  to  be  extremely  large 
(even  for  small  applied  displacements),  and  approximately  invariant  along 
the  boundary. 

The  large  magnitude  and  approximate  constancy  of  the  stress  at  the 

plastic  boundary,  and  the  similarity  between  the  boundary  trajectory  and 

the  observed  damage  zones  strongly  suggest  that  the  propagation  of  the 

damage  is  a consequence  of  the  onset  of  matrix  yield.  The  bundle  failure 

stress  a is  thus  likely  to  be  dictated  by  the  conditions  required  for 
c 

matrix  plasticity,  derived  from  the  slip  line  field  solutions 


a ” 0 / (h/a) 
c y 


where/  is  plotted  in  Fig.  4.  Combining  with  eqn.(l)  yields  the  critical 
applied  stress; 


a -o  /(h/a)  (E-  /E.  ^)/(l+d,  ^/d.  )‘ 

®c  y long  lat  lat  long 


The  extensive  plastic  flow  that  accompanies  damage  propagation  will 
cause  the  broken  fibers  to  displace  plastically  (Fig.  5a).  Hence,  when 
the  applied  stress  is  removed,  the  residual  elastic  strain  in  the  fibers 
coupled  with  the  constraint  at  the  damage  site  (Fig.  5b)  could  cause  a 
second  series  of  fiber  fractures.  This  type  of  fracture  is  fre>^uently 
observed  within  the  damage  zone  (Fig.  1). 
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The  observations  and  analysis  of  lateral  bundle  damage  contain  some 
specific  inferences  about  microstructural  modifications  that  would  enhance 
the  damage  resistance.  The  first  important  parameter,  (as  for  kinkingl) 
is  the  matrix  yield  strength, : the  critical  damage  stress  being  directly 
proportional  to  a^.  The  relative  dimensions  of  the  lateral  bundles  are 
also  important.  The  pressure  required  for  matrix  yield  increases  as  h/a 
increases  (through  f(h/a))  and  hence,  a small  ratio  or  the  lateral 
dimension  to  the  longitudinal  dimension  of  the  lateral  bundles  (up  to 
1/10)  would  be  of  significant  benefit.  The  importance  of  the  relative 
dimensions  of  the  lateral  and  longitudinal  bundles  depends  on  the  response 
of  the  latter.  Provided  that  the  restraint  of  the  longitudinal  bundles  is 
intact,  eqn.(8)  indicates  that  a large  proportion  of  longitudinal  bundles 
is  preferred.  Since  a large  longitudinal  bundle  content  also  leads  to  a 
greater  kink  resistance,!  this  is  a generally  desirable  microstructure. 
However,  it  shoud  be  recognized  that  once  kinking  of  the  longitudinal 
bundles  initiates  in  a material  with  a large  content  of  these  bundles, 
extensive  lateral  bundle  damage  would  quickly  ensue. 

The  effect  of  an  orthogonal  stress  applied  along  the  axis  of  the 
lateral  bundles  is  also  of  interest,  especially  for  impact  damage 
analysis. 2 in  terms  of  the  slip  line  field  solution  (Fig.  3),  the 
orthogonal  stress  would  supJSly  counteracting  shear  stress  along  45° 
trajectories.  However,  even  for  hydrostatic  compression,  shear  stress 
would  exist  outside  the  45°  zone  AQA* , and  damage  could  still 
occur.  A relatively  small  dependence  of  the  damage  condition  on  the 
orthogonal  stress  would  thus  be  predicted. 
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Estimation  of  Bundle  Displacement 

The  displacements  at  (r,9)  in  a semi-infinite  body  subjected  to  a line 
force  P per  unit  length  are; 


_ 2P 

ir<E> 


(i-nr/z*)  COS0  -f (1-v)  9sin0 
2 


■n’<E> 


(l+v)sin0  + (Jln(r/z*))Bin0  - 


(Al) 


where  r and  9 are  the  position  coordinates.  Expressions  (Al)  can  be  used 
to  obtain  the  displacement  due  to  a distribution  of  pressure,  a (x) . The 
horizontal  displacement,  u,  close  to  the  surface  (0'''IT/2)  are  relatively 
uniform  and  have  little  effect  on  the  stresses  in  the  fibers.  Hence, 
solutions  are  only  derived  for  the  vertical  displacements,  v. 

Consider  the  displacements  at  a point  (x,z)  due  to  a small  element  d a 
of  pressure.  The  vertical  displacement  from  eqn(Al)  is: 

dv(z,x)  =»  |^(l+v)sin0-(l-v)9cos9  + 2sin0  !ln(z/z*cos0)J  da 

(A2) 
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tan 


-1  f 


v(z,x) 


az 


ir<E>  , 

.iJ  ' 

-tan  (x/ z) 


[(l+v)  t an  2/2.]  slsf  . IbiJi  . 2sl«  Jncose  de 
COS^e  COS0  cos‘^0  I 


u<E> 


(1-V) 


p 

/ \ 

J 2,,,  ,2  . 

V z +(djj^-x)  + '^z  +x 

+ 2 '^z^+Cdj^-x)^  An  1 

^z*/v4^+(dj^-x)^j 

+ 2 >/  z^4x^  iln(z*/N4^+x^  )+  z(l- 


( ( -l/x\^^  / 

ytan  +ytan  [-)/ 

^ O 


.144/ 


Along  Che  surface,  z=0,  this  reduces  to: 


(dj^-x)  An 


z*/(dj^-x) 


+ x£,n(z*/x)  + dj^(l-v)/2| 


(A4) 
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FIGURE  CAPTIONS 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Polarized  ligiit  reflected  micrographs  of  complex  kinking 
damage  (a)  a typical  damage  zone  showing  the  trajectory 
of  the  damage  between  opposing  comers  of  the  out-of-plane 
bundles  (b)  an  atypical  damage  zone. 

A schematic  indicating  the  fiber  displacement  in  response 
to  the  local  stress  G 

Slip  line  fields  for  compression  of  a plate  between 
plattens  (a)  aligned  plattens  (b)  misaligned  plattens. 
The  indentation  pressure  as  a function  of  the  relative 
dimensions  h/a  of  the  lateral  bundles. 

A schematic  indicating  a probable  mode  of  double 
fracture:  (a)  the  kink  fracture  when  the  force  is 

applied,  and  the  permanent  plastic  displacement  at  the 
damage  zone;  (b)  the  elastic  relaxations  and  second 
fracture  that  occur  at  the  region  of  high  tension  when  the 
applied  force  is  removed. 
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